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Abstract 

In this paper, we introduce a new class of distributions which is obtained by compounding the extended 
Weibull and power series distributions. The compounding procedure follows the same set-up carried 
out by Adamidis and Loukas (1998) and defines at least new 68 sub-models. This class includes some 
well-known mixing distributions, such as the Weibull power series (Morais and Barreto-Souza, 2010) and 
exponential power series (Chahkandi and Ganjali, 2009) distributions. Some mathematical properties of 
the new class are studied including moments and generating function. We provide the density function of 
the order statistics and obtain their moments. The method of maximum likelihood is used for estimating 
the model parameters and an EM algorithm is proposed for computing the estimates. Special distribu- 
tions are investigated in some detail. An application to a real data set is given to show the flexibility 
and potentiality of the new class of distributions. 

Keywords: EM algorithm, Extended Weibull distribution, Extended Weibull power series distribution, 
Order statistic, Power series distribution. 



1. Introduction 

The modeling and analysis of lifetimes is an important aspect of statistical work in a wide variety 
of scientific and technological fields. Several distributions have been proposed in the literature to model 
lifetime data by compounding some useful lifetime distributions. Adamidis and Loukas (1998) introduced 
a two-parameter exponential-geometric (EG) distribution by compounding an exponential distribution 
with a geometric distribution. In the same way, the exponential Poisson (EP) and exponential logarith- 
mic (EL) distributions were introduced and studied by Kus (2007) and Tahmasbi and Rezaei (2008), 
respectively. Recently, Chahkandi and Ganjali (2009) proposed the exponential power series (EPS) fa- 
mily of distributions, which contains as special cases these distributions. Barreto-Souza et al. (2010) 
and Lu and Shi (2011) introduced the Weibull-geometric (WG) and Weibull-Poisson (WP) distributions 
which naturally extend the EG and EP distributions, respectively. In a very recent paper, Morais and 
Barreto-Souza (2011) defined the Weibull power series (WPS) class of distributions which contains the 
EPS distributions as sub-models. The WPS distributions can have an increasing, decreasing and upside 
down bathtub failure rate function. 

Now, consider the class of extended Weibull (EW) distributions, as proposed by Gurvich et al. (1997), 
having the cumulative distribution function (cdf) 

G(x; a,€) = l- e - aH( - x '<t\ x>0, a > 0, (1) 

where H (x; £) is a non-negative monotonically increasing function which depends on a parameter vector 
£. The corresponding probability density function (pdf) is given by 

g{x; a,£) =ah(x; £) e""^ 2 ^, x > 0, a > 0, (2) 

where h(x; £) is the derivative of H(x; £). 
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Note that many well-known models are special cases of equation ([I]) such as: 
(i) H(x;£) — x gives the exponential distribution; 

(tt) H(x;£) = x 2 yields the Rayleigh distribution (Burr type-X distribution); 

(Hi) H(x;£) = \og(x/k) leads to the Pareto distribution; 

(iv) H(x;£) — /3 _1 [exp(/3x) — 1] gives the Gompertz distribution. 

In this article, we define the extended Weibull power series (EWPS) class of univariate distributions 
obtained by compounding the extended Weibull and power series distributions. The compounding pro- 
cedure follows the key idea of Adamidis and Loukas (1998) or, more generally, by Chahkandi and Ganjali 
(2009) and Morais and Barreto-Souza et at (2011). The new class of distributions contains as special 
models the WPS distributions, which in turn extends the EPS distributions and defines at least new 68 
(17 x 4) sub-models as special cases. The hazard function of our class can be decreasing, increasing, 
bathtub and upside down bathtub. 

We are motivated to introduce the EWPS distributions because of the wide usage of the general class 
of Weibull distributions and the fact that the current generalization provides means of its continuous 
extension to still more complex situations. 

This paper is organized as follows. In Section 2, we define the EWPS class of distributions and 
demonstrate that there are many existing models which can be deduced as special cases of the proposed 
unified model. In Section 3, we provide the density, survival and hazard rate functions and derive some 
useful expansions. In Section 4, we obtain its quantiles, ordinary and incomplete moments. Further, 
the order statistics are discussed and their moments are determined. Section 5 deals with reliability and 
average lifetime. Estimation of the parameters by maximum likelihood using an EM algorithm and large 
sample inference are investigated in Section 6. In Section 7, we present suitable constraints leading to 
the maximum entropy characterization of the new class. Three special cases of the proposed class are 
studied in Section 8. In Section 9, we provide an application to a real data set. The paper is concluded 
in Section 10. 



2. The new class 

Our class can be derived as follows. Given N, let X\,... be independent and identically dis- 
tributed (iid) random variables following ([I]). Here, N is a discrete random variable following a power 
series distribution (truncated at zero) with probability mass function 

p n = P[N = n) = -g^-,n = l,2,..., (3) 

where a n depends only on n, 

c i e ) = ££Li a n n and 6> > is such that C{6) is finite. Table [l] 
summarizes some power series distributions (truncated at zero) defined according to ([3| such as the 
Poisson, logarithmic, geometric and binomial distributions. Let Xm = rnin{Xi}^ =1 . The conditional 
cumulative distribution of Xn\\N = n is given by 



G X(1) \N=n( x ) — 1 



C 



-naH (x;£) 



i.e., X(i)|iV = n follows a general class of distributions with parameters na and £ based on the same 
H(x;£) function. Hence, we obtain 



P(X {1) <x,N = n) 



x > 0, n > 1. 



C(0) 

The EWPS class of distributions can then be defined by the marginal cdf of Xm; 

F(x;9,a,i) = l- [ J , x > 0. (4) 
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Table 1: Useful quantities for some power series distributions. 

The random variable X following Q with parameters 9 and a and the vector £ of parameters is 
denoted by X ~ EWPS(#,a,£). Equation Q extends several distributions which have been studied in 
the literature. The EG distribution (Adamidis and Loukas, 1998) is obtained by taking H[x\ £) = x 
and C{9) = 0(1 - A)" 1 with t? G (0, 1). Further, for H(x; £) = x, we obtain the EP (Kus, 2007) and EL 
(Tahmasbi and Rezaei, 2008) distributions by taking C{9) = e e - 1, 6 > 0, and C{9) = - log(l -6),0e 
(0, 1), respectively. In the same way, for H(x; £) = x 7 , we obtain the WG (Barreto-Souza et ai, 2009) 
and WP (Lu and Shi, 2011) distributions. The EPS distributions are obtained from Q by mixing 
H{x\ £) = x with any C(9) listed in Table [l] (see Chahkandi and Ganjali, 2009). Finally, we obtain the 
WPS distributions from Q by compounding H(x; £) = x 1 with any C{9) in Table [T] (see Morais and 
Barreto-Souza, 2011). Table [2] displays some useful quantities and respective parameter vectors for each 
particular distribution. 

3. Density, survival and hazard functions 

The density function associated to Q is given by 

C" (ft p —aH(x; £)\ 

f{x; e,a,H) = 9a h(x; g) e~ aH ^ Q - 1 — >-, x > 0. (5) 

Proposition 1. The EW class of distributions with parameters col and £ is a limiting special case of 
the EWPS class of distributions when 9 — > + , where c — min{n £ N : a n > 0}. 

Proof. This proof uses a similar argument to that found in Morais and Barreto-Souza (2011). Define 
c = min {n G N : a n > 0}. We have 



lim F(x) = 1 — lim 

0^0+ S-KI+ 



oo 

OO 

£a„ 9 n 



e -caH(x;£) + a -l ^ ^ jn-Cg-naHfefl 
n— c+1 



= 1 — lim 

1 + a- 1 a n 9 n - c 

n— c+1 

for X > 0. □ 



We now provide an interesting expansion for ([5]). We have C'{9) = YlnLi n a « ^™ By using this 
result in ([5]), it follows that 

oo 

f(x;6,a,£) = ^2p n g(x; na,£), (6) 

71=1 

where c/(x; na, £) is given by ([2]). Based on equation ([6]), we obtain 

oo 

F(x;6,a,Z) = l-J2Pne- naH{x; V. 
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Hence, the EWPS density function is an infinite mixture of EW densities. So, some mathematical 
quantities (such as ordinary and incomplete moments, generating function and mean deviations) of 
the EWPS distributions can be obtained by knowing those quantities for the baseline density function 
g(x; na,£). 

The EWPS survival function is given by 



C (0 e -aJI(xi€)) 



and the corresponding hazard rate function becomes 



t(x; 9, a, £) = 9a h(x; £) e 



C{9) 



(7) 



(8) 



4. Quantiles, moments and order statistics 

The EWPS distributions are easily simulated from Q as follows: if U has a uniform C/(0, 1) distri- 
bution, then the solution of the nonlinear equation 



X = H~ 



log 



C-\C(9)(l-U)) 



has the EWPS(6>, a, £) distribution, where H~ 1 (-) and C _1 (-) are the inverse functions of H(-) and 
C(-), respectively. To simulate data from this nonlinear equation, we can use the matrix programming 
language Ox through SolveNLE subroutine (see Doornik, 2007). 

We now derive a general expression for the rth raw moment of X, which may be determined by using 
([6| and the monotone convergence theorem. So, for r G N, we obtain 



(9) 



where Z is a random variable with pdf g(z; na, £). 

The incomplete moments and moment generating function (mgf ) follow by using ^ and the monotone 
convergence theorem: 



ry 00 

lx(y) = x r f{x)dx = ^2p n I z (y) 

J n=l 



and 



where Z is defined as before. 

Order statistics are among the most fundamental tools in non-parametric statistics and inference. 
They enter in the problems of estimation and hypothesis tests in a variety of ways. Therefore, we now 
discuss some properties of the order statistics for the proposed class of distributions. The pdf fi :m (x) of 
the ith order statistic for a random sample Xi, . . . , X m from the EWPS distribution is given by 



fi:m{x*) 



mi 



(i- l)!(m-t)! 



f{x;0,a,£) 



' C{9) 



C{6) 



> 0, (10) 



where /(x; 9, a, £) is the pdf given by By using the binomial expansion, we can write ( 10 ) as 

■ i— 1 /■ -,\ 



fi:m{x) , . 



(i - l)\(m-i)\ 



f(x;9,a,£)J2(-l) 



3=0 
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(ii) 
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where S(x; 0, a, £) is given by Q. The corresponding cumulative function is 



oo m /U\ / 



An alternative form for (10) can be obtained from ^ as 



fi:m (#) 



(i- l)!(m-i)! ^ 

K J v y n— 1 j— 



]T ^ Uj Pn g(x; na, t)S(x; 9, a, 



(12) 



where Wj = (* 1 ). So, the sth raw moment Xi :m comes immediately from the above equation 



oo i—1 



(13) 



(i- l)!(m-i)! 

v ; v ; n=l j=0 

where Z ~ EW(na,^) is defined before. 

5. Reliability and average lifetime 

In the context of reliability, the stress-strength model describes the life of a component which has a 
random strength X subjected to a random stress Y. The component fails at the instant that the stress 
applied to it exceeds the strength, and the component will function satisfactorily whenever X > Y. 
Hence, R — P(X > Y) is a measure of component reliability. It has many applications, especially in 
engineering concepts. The algebraic form for R has been worked out for the majority of the well-known 
distributions. Here, we obtain the form for the reliability R when X and Y are independent random 
variables having the same EWPS distribution. 

The quantity R can be expressed as 



/*OG 

R= I f(x;9,a,£)F(x;6,a,€)dx. 
Jo 

Substituting Q and ^ into equation (14), we obtain 

C'(fe- Qff <^)) 



(14) 



R 
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where the integral can be calculated from the baseline EW distribution. 
The average lifetime is given by 



oo 

oo 

5> J dx. 



Given that there was no failure prior to Xo, the residual life is the period from time xq until the time 
of failure. The mean residual lifetime can be expressed as 

oo 

m(x ; 9, a, £) = [Pi(X > xq)}^ 1 J y f(x + y; 9, a, £)dy 

o 

oo °° 

= [S^o^^Pn / yg{x + y;na,£)dy. 

n=l { 

The last integral can be computed from the baseline EW distribution. Furthermore, m(xo; 9, a, £) — > 
E(X) as x -> 0. 







6. Maximum likelihood estimation 



6. 1 . Preliminaries 



Here, we determine the maximum likelihood estimates (MLEs) of the parameters of the EWPS class 
of distributions from complete samples only. Let X\ , . . . , X n be a random sample with observed values 
xi,...,x n from an EWPS distribution with parameters 9, a and Let 9 = (9, a,£) T be the pxl 
parameter vector. The total log-likelihood function is given by 



= £ n (x; 9) = n [log 9 + log a - log C{9)] - a £ H{ Xi ; $) + ^ log h( Xi ; £) 

n 

+ ^logC'^e-"^ 5 ^). 



(15) 



The log-likelihood can be maximized either directly by using the SAS (PROC NLMIXED) or the Ox 
program (sub-routine MaxBFGS) (see Doornik, 2007) or by solving the nonlinear likelihood equations ob- 
tained by differentiating ( 15 1. The components of the score function U n (&) = (d£ n /d9, d£ n /da, d£ n /dt;) 
are 
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For interval estimation on the model parameters, we require the observed information matrix 



V Uee u atL 



} 



whose elements are listed in Appendix A Let O be the MLE of 8. Under standard regular conditions 
stated in Cox and Hinkley (1974) that are fulfilled for our model whenever the parameters are in the 

interior of the parameter space, we have that the asymptotic distribution of ^fn ^9 — 9^j is multivariate 

normal N p (0, K(Q)^ 1 ), where K(Q) — lim„_ i . 00 J n (Q) is the unit information matrix and p is the number 
of parameters of the compounded distribution. 



6.2. The EM algorithm 

Here, we propose an EM algorithm (Dempster et ai, 1977) to estimate 9. The EM algorithm is a 
recurrent method such that each step consists of an estimate of the expected value of a hypothetical 
random variable and then maximizes the log-likelihood for the complete data. Let the complete-data 
be X\, . . . , X n with observed values Xx,...,x n and the hypothetical random variables Z\, . . . ,Z n . The 
joint probability function is such that the marginal density of X%, . . . ,X n is the likelihood of interest. 
Then, we define a hypothetical complete-data distribution for each (Xi 1 Zi) T , i = 1, . . . ,n, with a joint 
probability function in the form 

9( X ,z;0)= a ^h(x; e)e-«*<*'«, 



7 



where 9 and a are positive, x > and z € N. Under this formulation, the E-step of an EM cycle requires 
the expectation of Z\X; 9< r ) = (flW^W^W) 1 as the current estimate (in the rth iteration) of 0. The 
probability function of Z given X, say <7(z|a;), is given by 

„( 7 \t\ - Za ^ 6 ~ 1 - a {z-l)H( Xi ;i) 

9[ZlX) - C'{6e- aH i*«t)) 

and its expected value is 

E{Z\X) = 1 + 9e~ aH ^ 6 [ - „ t J . 

The EM cycle is completed with the M-step by using the maximum likelihood estimation over 6, 
where the missing Z's are replaced by their conditional expectations given before. The log-likelihood for 
the complete-data is 



t n {x u ...,x n ; Z!,.. .,z n ; a, 6, £) cx nlog a + log 6 ^ z, + ^ log h(x,; £) 

»=1 i=l 

n 

-aj]^(ii; £)-nlogC(0). 

i=l 

So, the components of the score function £/*(0) = (<9Z* 9Z* /9a, 9i* /<9£) T are 
<9Z* n C"(0) dl* n ^ rT . 



i=l v ' i=l 

and 

az* _ " aiog/,^;!) " agfo; g) 

% ^ ^ a h Zl d t* ' 

From a nonlinear system of equations U*(Q) = 0, we obtain the iterative procedure of the EM 
algorithm 

a <«-i) = ? = ^ (t+1) ) i v 



and 



^ aio g fe(si;ff t+1) ) , (t) ^ (p gg^gw)) _ 

where d( t+1 ) and |( t+1 ) are obtained numerically. Here, for i = 1, . . . , n, we have 

~, s n» ( a(t) „-<i (t) Hfxn £ (t) )\ 

Note that, in each step, 9, a and £ are estimated independently. The EWPS distributions can be 
very useful in modeling lifetime data and practitioners may be interested in fitting one of our models. 



7. Maximum entropy identification 

Shannon (1948) introduced the probabilistic definition of entropy which is closely connected with the 
definition of entropy in statistical mechanics. Let X be a random variable of a continuous distribution 
with density /. Then, the Shannon entropy of X is defined by 

Wsh(f) = ~ I f(x; 9, a, £) log [f(x; 9, a, £)] dx. (16) 
Jr 
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Jaynes (1957) introduced one of the most powerful techniques employed in the field of probability and 
statistics called the maximum entropy method. This method is closely related to the Shannon entropy 
and considers a class of density functions 



F = {f(x; 9, a, £) : E f (T t (X)) = a u i = 0, . . . , m} , (17) 

where Ti(X),i = 1, ...,m, are absolutely integrable functions with respect to /, and T (X) = a = 1. 
In the continuous case, the maximum entropy principle suggests deriving the unknown density function 
of the random variable X by the model that maximizes the Shannon entropy in (16), subject to the 
information constraints defined in the class F. Shore and Johnson (1980) treated axiomatically the 
maximum entropy method. This method has been successfully applied in a wide variety of fields and has 
also been used for the characterization of several standard probability distributions; see, for example, 
Kapur (1989), Soofi (2000) and Zografos and Balakrishnan (2009). 

The maximum entropy distribution is the density of the class F, denoted by f ME , which is obtained 
as the solution of the optimization problem 

f ME {x;9,a,£) = argmaxHsh. 

Jaynes (1957, p. 623) states that the maximum entropy distribution f , obtained by the constrained 
maximization problem described above, "is the only unbiased assignment we can make; to use any other 
would amount to arbitrary assumption of information which by hypothesis we do not have". It is the 
distribution which should not incorporate additional exterior information other than which is specified 
by the constraints. 

We now derive suitable constraints in order to provide a maximum entropy characterization for our 
class of distributions defined by ([4]). For this purpose, the next result plays an important role. 

Proposition 2. Let X be a random variable with pdf given by Then, 

CI. E[log{C'(9e- aH( - x -^))] = -^-E[C'(0e- aff ( F; «))log(C"(0e- Qff ( y '«'))l ; 

C{6) 

C2. Epog(ft(*; €))] = ^-E[C'(ee- H ^)log(h(Y; £))] ; 

C3. E [H(X; Q] = JLe [C '(9 e - aH ( Y ^)H(Y ; £)] , 

where Y follows the EW distribution with density function pj). 

Proof. The constraints CI, C2 and C3 are easily obtained and therefore their demonstrations are omitted. 

□ 

The next proposition reveals that the EWPS distribution has maximum entropy in the class of all 
probability distributions specified by the constraints stated in the previous proposition. 

Proposition 3. The pdf f of a random variable X, given by ([5]), is the unique solution of the optimization 
problem 

f(x;6,a,£) = argmaxHsTj, 

h 

under the constraints CI, C2 and C3 presented in the Proposition^ 

Proof. Let r be a pdf which satisfies the constraints CI, C2 and C3. The Kullback-Leibler divergence 
between r and / is 

D(r,f)= [ T(x;6,a,®\og( T p d ' a 'f )dx. 
Jr \f[x;9,a,£) J 

Following Cover and Thomas (1991), we obtain 

< D(t, f) = [ t(x; 9, a, £) log [r(x; 9, a, £)] dx- t(x; 9, a, £) log [f(x; 9, a, £)] dx 
Jr Jr 

= -m Sh (T]0,a,t)- ( T(x;e,a,£)\og[f{x;6,a,t)]dx. 
Jr 



From the definition of / and based on the constraints CI, C2 and C3, it follows that 



t(x) log [f(x)j dx = log(6a) 



C(6) 



e{c'(9 e- aH ( Y -> O ) log [h(Y; £)] } - log [C{6)] 



a 



C(9) 



E 



C"(9e- aH(Y <^)H(Y; £ 



+ ^yE{log [C(ee- aH ^)] C'(9e- aH ™)} 

= [ f(x; 9, a, log [f{x- 9, a, £)] dx = -M sh (f), 
Jw 



where Y is defined as before. So, we have Hs/^t) < W-Sh(f) with equality if and only if r(x;9,a,^) = 
f(x; 9, a, £) for all x, except for a set of measure 0, thus proving the uniqueness. □ 

The intermediate steps in the above proof in fact provide the following explicit expression for the 
Shannon entropy of the EWPS distribution 



!W/) = -log(0a) 
9 

E 



C(8) 



C{8) 

C'(9e~ aH{Y ^ ) )H(Y; i) 



E {C"(0e- Q "( y >«) log [h(Y; £)]} + log [C{6)\ 

E [C{e e ) log [C (9 e~« H ^ «))]}. (18) 



For some EWPS distributions, the above results can only be obtained numerically. 



8. Special models 

In this section, we investigate some special cases of the EWPS class of distributions. We offer 
some expressions for moments and moments of the order statistics. To illustrate the flexibility of these 
distributions, we provide plots of the density and hazard rate functions for selected parameter values. 

8.1. Modified Weibull geometric distribution 

The modified Weibull geometric (MWG) distribution is defined by the cdf Q with H(x; £) = x 1 and 
C(9) = 6(1 - 9)- 1 leading to 

(1 - 9) exp (-ra^e Al ) , v 

^;M,7,A) = l- \_; exp P ( ^ 7eA ^ , *>0, (19) 

where 9 £ (0, 1). The associated pdf and hazard rate function are 

_ t exp (Xx — ax 7 e Ax ) 



and 



f(x;8,a, 1 ,\)=a(l-6){ 1 + \x)xi- 

[1 - 6>cxp (~ax^e Xx )} 



-(x; 9, a, 7 , A) = a( 7 + Xx) x^ 1 e M^L 



1 - 9 exp (-ax~'e Xx ) 

for x > 0, respectively. The MWG distribution contains the WG distribution (Barreto-Souza et al. 
(2010)) as the particular choice A = 0. Further, for A = and a = 1, we obtain the EG distribution 
(Adamidis and Loukas (1998)). Figures [l] and [2] display the density and hazard functions of the MWG 
distribution for selected parameter values. 

The rth raw moment of the random variable X having the MWG distribution has closed-form. It is 
calculated from ^ as 

oo 

E(X r ) = Y,Pn^r(n), (20) 

71=1 
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(a) a = 2, 7 = 1.5 and A = 0.5 



(b) a = 2, 7 = 0.8 and A = 0.01 



(c) a = 0.1, 7 = 6 and A = 0.5 



Figure 1: Plots of the MWG density functions for 9 = 0.01 (solid line), 9 = 0.2 (dashed line), 9 = 0.5 (dotted line) and 
9 = 0.9 (dotdash line). 



where [i r (n) = J °° x r g(x; na, 7, X)dx denotes the rth raw moment of the MW distribution with parame- 
ters na, 7 and A. Here p n corresponds to the probability function of the geometric distribution. Carrasco 
et al. (2008) determined an infinite representation for the rth raw moment of the MW distribution with 
these parameters expressed as 



MrO) = 

i\ ,. . . ,l r — 1 



4i 1 ,...,i r r(a r /7 + 1) 
(na) s rh 



(21) 



where 
and 



■Ail,...,i r — Oil ? ■ • • i a i r a nd S r — i\, . . . , i r 



(-1)*+V- 2 fX 



(i-iy. Vt 



Hence, the moments of the MWG distribution can be obtained directly from equations (20) and (21 1 




0.0 0.5 1.0 1.5 0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 

XXX 

(a) a = 2, 7 = 1.5 and A = 0.5 (b) a = 2, 7 = 0.8 and A = 0.01 (c) a = 0.1, 7 = 6 and A = 0.5 

Figure 2: Plots of the MWG hazard rate functions for 9 = 0.01 (solid line), 9 = 0.2 (dashed line), 9 = 0.5 (dotted line) and 
9 = 0.9 (dotdash line). 

The density of the ith order statistic X i:m in a random sample of size m from the MWG distribution 
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is given by (for 1 = 1,..., m) 



fi:m(x) = T . T~\T7~~ ^7 Y] Y] OJj p n 

(i — 1 ) (m — in 

y ' v ' n=l 7=0 



(1 — 0) cxp (— ax^e 
1 - 6exp (-axie Xx ) 



- m-\-j — i 



g(x;na,j,X), 



where g(x; na, 7, A) denotes the MW density function with parameters na, 7 and A. From ( 13 ), we obtain 



(i-1) (m-i) jyn ^ 



(l-6») cxp (-a jr>c AX ) 
1 - 6»exp(-aXTe AX ) 



- m+j — i 



8.2. Pareto Poisson distribution 



The Pareto Poisson (PP) distribution is defined by taking H(x; £) = log(x/fc) and C{ff) = e e — 1 
in Q, which yields 

„ ,n exp \$ (k/x) a ] — 1 
F{x;9,a,k) = 1 V 1 > ^ ^ fc - 



e e - 1 



The pdf and hazard functions of the PP distribution are 

6ak a exp[9{k/x) a ] 



and 



f(x;6,a,k) 
t(x; 0, a, k) 



<9a/fc Q exp [(9(/c/:r) Q ] 



a^+Hexp [0 (fc/a;) a ] - 1}' 

We obtain the Pareto distribution as a sub-model when — > 0. The rth moment of the random 
variable X following the PP distribution becomes 



E{X r ) 



(e e - 1) ^ (n - 1)! (na - r) ' 



na > r. 



(22) 



In particular, setting r = 1 in (22|, the mean of X reduces to 

ik 



an v - 

- 1 ^ (n- 1 



' (n - 1)! (na - 1) 



, na > 1. 
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(a) a — k — 0.5 



(b) a = 2 and k = 1 



(c) a = 7 and k = 2 



Figure 4: Plots of the PP hazard functions for 9 = 0.01 (solid line), 9 = 0.2 (dashed line), 9 = 0.5 (dotted line) and 9 = 0.9 
(dotdash line). 



From equation (13 1, the sth moment of the ith order statistic, for i = 1, . .. ,m, is given by 



E 



( ^™ } = (i-l)7(m-i)l £ £ u i E 
v ' v ' n=l 7=0 



exp(fl(fc/X) a )-l 
e e - 1 



m+j— i 



where p n denotes the Poisson probability function. Furthermore, after some algebra, the Shannon entropy 
for the PP distribution reduces to 



BW/) = log 



where 
//I = E exp ^9 

fi2 = E exp 

and 

Ms = E 
where 



log 
log 



c g - 1 
9a 



1 



c e - 1 



(pi - + M3) • 



Chi(20) - log(20) + Shi(20) - 7 



2(e e - 1) L " 
Chi(20) - log(2fl) + Shi(2fl) - 7 
2a(e e - 1) 



(e^-l)logfe , 



a0fc 



4(e 8 - 1) 



!l{l-(20 + l)e 29 } 



< cosh(i) — 1 
Chi(z) = 7 + log z + I — dt 



is the hyperbolic cosine integral, 



Shi(z) 



sinh(i) - 1 



dt 



is the hyperbolic sine integral and 7 ~ 0.577216 is the Euler-Mascheroni constant. 

8.3. Chen logarithmic distribution 

The Chen logarithmic (CL) distribution is defined by the cdf Q with H(x; £) = cxp(x^) — 1 and 
C{9) = - log(l - 9), leading to 

FM _, log{l-eexp[-a(exp(^)-l)]} 
W_ log(l-0) ' ' 

13 



where 9 £ (0, 1). The associated pdf and hazard rate function (for x > 0) are 

0a.bx b - 1 exp {x b - a [exp(a; b ) - l] } 
" X ' = log(l - 9) {6 exp [-a(exp(x b ) - 1)] - 1} 

and 

. Qabx b ~ x exp [a; 6 - a(exp(x 6 ) - 1)1 

t(x) = 



{6»exp[-a(exp(a; fc ) - 1)] - l}log{l - 6 exp [-a(exp(a; b ) - 1)]}' 
respectively. 





~i 1 1 r 

0.4 0.6 0.8 1.0 





(a) a = 2 and 6=1 



(b) a = 3.5 and 6 = 7 



(c) a = 3.5 and 6 = 1.5 



Figure 6: Plots of the CL hazard rate functions for 8 = 0.01 (solid line), 8 = 0.2 (dashed line), 8 = 0.5 (dotted line) and 
8 = 0.9 (dotdash line). 

The density of the ith order statistic X^ m in a random sample of size m from the CL distribution is 
given by (for i = 1, . . . , m) 



OO I— 1 

uj* p n g(x; na, b) |log 1 — exp(a — a e x ) j 



m-\-j — 1 



u — 1 ) ( m — z ) 

v ' K ' n=l j=0 

where g(x; na, b) is the pdf of the Chen distribution with parameters na and b and p n denotes the 
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logarithmic probability mass function and 



cr) 


i 






Liog(i-0)J 



In the same way, the sth raw moment of X^ :m is obtained directly from 

| oo i-1 

E (X? :m ) = -— ^— £ Pn E i ZS ex P [Mm +3 ~ 1)(1 - exp(Z b ))] } , 
^ n=l j=0 

where Z ~ Chen(na, 6). 
9. Application 

Fonseca and Franga (2007) studied the soil fertility influence and the characterization of the biologic 
fixation of N2 for the Dimorphandra wilsonii rizz growth. For 128 plants, they made measures of the 
phosphorus concentration in the leaves. The data are listed in Table [3j We fit the MWG, Gompertz 
Poisson (GP), PP, Chen Poisson (CP) and CL models to these data. We also fit the three-parameter 
WG distribution introduced by Barreto-Souza et al. (2010). The required numerical evaluations are 
implemented using the SAS (PROCNLMIXED) and R softwares. 



0.22 


0.17 


0.11 


0.10 


0.15 


0.06 


0.05 


0.07 


0.12 


0.09 


0.23 


0.25 


0.23 


0.24 


0.20 


0.08 


0.11 


0.12 


0.10 


0.06 


0.20 


0.17 


0.20 


0.11 


0.16 


0.09 


0.10 


0.12 


0.12 


0.10 


0.09 


0.17 


0.19 


0.21 


0.18 


0.26 


0.19 


0.17 


0.18 


0.20 


0.24 


0.19 


0.21 


0.22 


0.17 


0.08 


0.08 


0.06 


0.09 


0.22 


0.23 


0.22 


0.19 


0.27 


0.16 


0.28 


0.11 


0.10 


0.20 


0.12 


0.15 


0.08 


0.12 


0.09 


0.14 


0.07 


0.09 


0.05 


0.06 


0.11 


0.16 


0.20 


0.25 


0.16 


0.13 


0.11 


0.11 


0.11 


0.08 


0.22 


0.11 


0.13 


0.12 


0.15 


0.12 


0.11 


0.11 


0.15 


0.10 


0.15 


0.17 


0.14 


0.12 


0.18 


0.14 


0.18 


0.13 


0.12 


0.14 


0.09 


0.10 


0.13 


0.09 


0.11 


0.11 


0.14 


0.07 


0.07 


0.19 


0.17 


0.18 


0.16 


0.19 


0.15 


0.07 


0.09 


0.17 


0.10 


0.08 


0.15 


0.21 


0.16 


0.08 


0.10 


0.06 


0.08 


0.12 


0.13 







Table 3: Phosphorus concentration in leaves data set. 



Min. Qi Q 2 Mean Q 3 Max. Var. 

0.0500 0.1000 0.1300 0.1408 0.1800 0.2800 0.0030 



Table 4: Descriptive statistics. 

Tables |4] and [5] display some descriptive statistics and the MLEs (with corresponding standard errors 
in parentheses) of the model parameters. Since the values of the Akaike information criterion (AIC), 
Bayesian information criterion (BIC) and consistent Akaike information criterion (CAIC) are smaller for 
the CL distribution compared with those values of the other models, this new distribution seems to be 
a very competitive model for these data. 

Plots of the estimated pdf and cdf of the MWG, WG, GP, PP, CP and CL models fitted to these data 
are displayed in Figure [7| They indicate that the CL distribution is superior to the other distributions 
in terms of model fitting. 

Table [(3] lists the values of the Kolmogorov-Smirnov (K-S) statistic and the values of —2£(Q). From 
these figures, we conclude that the CL distribution provides a better fit to these data than the MWG, 
WG, GP, PP and CP models. 

10. Concluding remarks 

We define a new lifetime class of distributions, called the extended Weibull power series (EWPS), 
which generalizes the Weibull power series class of distributions proposed by Morais and Barreto-Souza 
(2011), which in turn extends the exponential power series class of distributions (Chahkandi and Ganjali, 
2009). We provide a mathematical treatment of the new distribution including expansions for the density 
function, moments, generating function and incomplete moments. Further, explicit expressions for the 
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(0.0017) 
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CP 
CL 


15.4386 
(22.8318) 
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(0.0001) 


14.7817 
(28.1576) 
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(0.0000) 


2.9212 
(0.2634) 
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Table 5: MLEs of the model parameters, the corresponding SEs (given in parentheses) and the statistics AIC, BIC and 
AICC. 
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Figure 7: Estimated (a) pdf and (b) cdf for the CL, MWG, PP, WG, GP and CP models to the percentage of Phosphorus 
concentration in leaves data. 



1G 



Model 


K-S 


-2£(6) 


MWG 


0.0943 


-393.6 


WG 


0.0873 


-384.5 


Gompcrtz Poisson 


0.1201 


-374.7 


Pareto Poisson 


0.3017 


-374.7 


Chen Poisson 


0.1159 


-389.7 


Chen Logarithmic 


0.0678 


-401.8 



Table 6: K-S statistics and — 2£(0) for the exccedances of phosphorus concentration in leaves data set. 



order statistics and Shannon entropy are derived. The EWPS density function can be expressed as 
a mixture of EW density functions. This property is important to obtain several other results. Our 
formulas related with the EWPS model are manageable, and with the use of modern computer resources 
with analytic and numerical capabilities, they may turn into adequate tools comprising the arsenal of 
applied statisticians. The estimation of the model parameters is approached by the method of maximum 
likelihood using the EM algorithm. The observed information matrix is derived. Further, maximum 
entropy identification for the EWPS distributions was discussed and some special models arc studied in 
some detail. Finally, we fit the EWPS model to a real data set to show the usefulness of the proposed 
class. We hope that this generalization may attract wider applications in the literature of the fatigue 
life distributions. 
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Appendix A. 

The elements of the p x p information matrix J n (0) are 



n 



n 



C"{0) ( C'(0) Y 
C{6) {C(6)J 



'E 

»=i 



Z2i_ 
Zli 



H( Xl ; i)e 



-2aH{xi;i) 



1=1 

n 



n n , 2 \ 
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where z u = C'(9e- aH ^^),z 2i = C"{6c- aH ^^) and z 3t = C" {ee~ aH ^^), for i = 1, . . . , n. 
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